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We show that two, non interacting 2-level systems, immersed in a common bath, can become 
mutually entangled when evolving according to a Markovian, completely positive reduced dynamics. 



The role of quantum entanglement is of primary impor- 
tance in quantum information and computation theory. 
In recent years, a lot of research has been devoted to 
study how to entangle two systems by means of a direct 
interaction between them (see for instance 0, S 0, S 01 ) • 
In such a context, the presence of an environment, e.g. 
a generic noisy reservoir or a heat bath, is commonly 
thought as counteracting entanglement creation, because 
of its decohering and mixing-enhancing effects. 

However, a heat bath can also provide an indirect inter- 
action between otherwise totally decoupled subsystems 
and thus a means to entangle them. Indeed, this has been 
explicitly shown in a simple, exactly solvable model 0. 
There, correlations between two subsystems are estab- 
lished during a transient phase where the reduced dy- 
namics of the subsystems contains memory effects. 

Instead, in this paper, we study the possibility that en- 
tanglement be created by the bath during the Markovian 
regime through a purely noisy mechanism. We consider 
two, non interacting 2-level systems, weakly coupled to a 
common heat bath. We then start with a total Hamilto- 
nian of the form 
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where , H and Hb drive the dynamics of the two 
subsystems and the bath in absence of each other; the 
interaction term couples each subsystem independently 
with the bath, and can be taken of the form: 
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where ci, 02, 03 are the Pauli matrices. Notice that we 
allow the subsystems to interact with the bath through 
different operators V a , while any direct coupling among 
themselves has been excluded. 

In the weak-coupling limit 0, S H E3, El E3> the 
reduced dynamics of the two 2-level systems takes on 
a Markovian form. Assuming a factorized initial state 
p <8> pb, where p is a state of the two subsystems and 
Pb is an equilibrium state of the bath, p evolves in time 
according to a quantum dynamical semigroup of com- 
pletely positive maps with generator of the Kossakowski- 
Lindblad form: 



d t p(t) = -i[H,p(t)] + L[p(t)} 
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The unitary term is the commutator with an effective 
Hamiltonian, H = ffW + H ^ + , consisting of 
single system pieces, including bath induced Lamb shifts, 
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plus, possibly, a bath generated two-system coupling 
term 
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The dissipative contribution L[p(t)] is as follows, 
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with T a = (To. C3> 1 for a = 1,2,3, T a = 1 ® cr Q -3 for 
a = 4,5,6, and T> — T>' a positive 6x6 matrix which 
guarantees the complete positivity of the evolution. By 
writing 
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with 3x3 matrices A = A\ C = C' and B, L[p] assumes 
a form more amenable to a physical interpretation: 
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A generator of this form has been applied in quantum 
optics to describe the phenomenon of collective resonance 
fluorescence (e.g. see 01)- I n the above expression, the 
first two contributions are dissipative terms that affect 
the first, respectively the second, system in absence of the 
other. On the contrary, the last two pieces represent the 
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way in which the noise may correlate the two subsystems; 
this effect is present only if the matrix B is different from 
zero. 

Remark 1 From the rigorous derivation of Markovian 
master equations 0, @ i one knows that the hamiltonian 
terms @, (0 and the entries of the matrix T> in © 
contain integrals of two point time-correlation functions 
of bath operators: Tr[pB V a Vp(t)]. In particular, the 
matrices [H^ 2 ^] in (JSJ) and [Bij] in (JSJ do not vanish 
only if the bath state ps correlate bath-operators V a 
pertaining to different subsystems, that is, if the expec- 
tations Tt[pb V a Vg(t)] are nonzero when 1 < a < 3 
and 4 < (3 < 6. Only in this case, entanglement has 
a chance to be created by the action of the bath. Indeed, 
if If ( 12 ) = and B = 0, the two subsystems evolve in- 
dependently and initially separable states may become 
more mixed, but certainly not entangled. 

In order to check whether the reduced two-system den- 
sity matrix p gets entangled at time t because of the 
time-evolution generated by equation J3J), one can use 
the partial transposition criterion 0, 

if p(t) acted upon with the partial transposi- 
tion with respect to one of the two subsystems 
has negative eigenvalues, then it is entangled; 
in the A- dimensional case we are studying, 
also the reciprocal is true, namely if p(t) is 
entangled, then, partial transposition makes 
negative eigenvalues appear. 

In physical terms, the bath is not able to create entan- 
glement if and only if the partial transposition preserves 
the positivity of the state p(t) for all times. 

Remark 2 Strictly speaking, this criterion allows us 
to study the possibility of creating entanglement start- 
ing from separable initial states. When the initial state 
is already entangled, the partial transposition criterion 
cannot settle the question; in such cases, the analysis of 
the entangling power of the bath can only be addressed 
through the study of how entanglement measures evolve 
in time under dissipative reduced dynamics. This prob- 
lem requires a separate treatment and will not be ad- 
dressed here. 

We therefore take separable states as initial states: as 
we shall see, this is not really a limitation for the pur- 
pose of discussing the possibility of bath-induced entan- 
glement creation. Further, we can restrict our study to 
pure states; indeed, if the bath cannot create entangle- 
ment out of these, it will certainly not entangle their 
mixtures. In view of this, we will consider initial states 
of the form 



where {|ai)}, i = 1,2, are orthonormal bases in 

the two dimensional Hilbert spaces of the two subsys- 
tems. For sake of defmiteness, we will operate the partial 
transposition over the second factor with respect to the 
basis {|&i), |6 2 )}- 

One can act with the partial transposition on both 
sides of equation and recast the result as 



d t p(t) = -i H,p(t) + L[p(t)} 
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here, p(t) denotes the partially transposed matrix p(t), 
while H is a new Hamiltonian to which both the unitary 
and the dissipative term in © contribute 
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where S is the diagonal 3x3 matrix given by: S — 
diag(— 1, 1, — 1). The additional piece L[-] is of the 
form ©, but with a new matrix T> — ► S ■ T> ■ S , where 
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p(0) = |ai)<oi|(g)|6i)<6i| 



and T denotes full transposition, while H^ 12 ^ is the coef- 
ficient matrix in J5|. 

Remark 3 Although p(t) evolves according to a mas- 
ter equation formally of Kossakowski-Lindblad form, the 
new coefficient matrix T> need not be positive. As a con- 
sequence, the time-evolution generated by (|10f> may re- 
sult to be neither completely positive, nor positive and 
therefore may not preserve the positivity of the initial 
state p(0) = p(0). 

Notice that both the hamiltonian and the dissipative 
terms of the original master equation @ contribute to 
the piece L[- ] in 1|1U|) . the only one that can produce neg- 
ative eigenvalues. In particular, this makes more trans- 
parent the physical mechanism according to which a di- 
rect hamiltonian coupling ij( 12 ) among the two systems 
can induce entanglement: on p(t), H^ 12 ^ "acts" as a dis- 
sipative contribution, which in general does not preserve 
positivity. The entanglement power of purely hamilto- 
nian couplings have been extensively studied in the re- 
cent literature Q, 0, El 0, • Instead, in the following we 
shall concentrate our attention on whether entanglement 
can be produced by the purely dissipative action of the 
heat bath; henceforth, we shall ignore the contribution 
of the matrix H^ 12 ^ in T>. In other words, we shall take 
into account only baths for which the induced two-system 



(9) hamiltonian coupling in (jSJ) is vanishingly small. 
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Remark 4 If D is positive, then, the time-evolution 
generated by iJTUjl is completely positive; therefore, p(t) 
is positive at all times and entanglement is not created. 
Instances of baths for which this happens can easily be 
provided: 

1. B = 0: in such a case, V is positive since such are A 
and C T , due to the positivity of T>\ this corresponds 
to a bath that does not dynamically correlate the 
two subsystems; 

2. IZe(B) — 0: as before, I? is block-diagonal and thus 
positive; 

3. Im(Bl = and C T = C or A T = A: in the first 
case, V = T>, while in the second T> = T> T ; 

4. A T = A and C T = C: in this case V = (V+V T )/2. 

In the last three cases, despite the fact that the two sub- 
systems are now dynamically correlated by the bath, the 
effect is not sufficient for entanglement production. Fur- 
ther, notice that entanglement is not created also in baths 
for which the corresponding coefficient matrix T> can be 
written as a convex combination of matrices satisfying 
the previous conditions. 

In order to check the presence of negative eigenvalues 
in p(t), instead of examining the full equation (|10J) we 
find convenient to study the quantity 

£(*) = MP(*)M » ( 14 ) 

where ip is any 4-dimensional vector. Assume that an 
initial separable state p has indeed developed a negative 
eigenvalue at time t, but not before. Then, there exists 
a vector state and a time t* < t such that £(t*) — 
0, £{t) > for t < t* and £(t) < for t > t*. The 
sign of entanglement creation may thus be given by a 
negative first derivative of £(i) at t = t*. Moreover, by 
assumption, the state p(t*) is separable. Without loss of 
generality, one can set t* = and, as already remarked, 
restrict the attention to factorized pure initial states. 

In other words, the two subsystems, initially prepared 
in a state p(0) = p(0) as in © , will become entangled by 
the noisy dynamics induced by their independent inter- 
action with the bath if: 1) £{0) = and 2) d t £(0) < 0, 
for a suitable vector \ip), 

2 

iv>> = ^ i a *> • ( 15 ) 

i J = l 

Given @, condition 1) readily implies: ipn = 0. 

Remark 5 Note that entanglement creation can not be 
detected by looking at the sign of the first derivative of 
£ (t) unless the test vector \ijj) is entangled itself. Indeed, 
£ (t) is never negative for a separable \ip). Thus, both 



components ■012 and ipzi in <|15|l have to be different from 
zero, since otherwise \ip) becomes separable. 

Remark 6 When dt£(0) > for all choices of the ini- 
tial state p(0) and probe vector \ip), the bath is not able 
to entangle the two systems, since p remains positive. 
The treatment of the case dt£(0) = requires special 
care: in order to check entanglement creation, higher or- 
der derivatives of £, possibly with a time dependent \ip) 1 
need to be examined. 

In order to prove that indeed there are baths for which 
£ (0) = and dt£(0) is negative, let us first make the 
choice |ai) = |&i) = |+) and |a 2 ) = |6 2 ) = |— ), where |±) 
are the eigenstates of 173; the general case is considered 
below. For = (|+) ® |-) + |-) ® |+>)/\/2, one finds 

d t £(0) = Tr[VU] , (16) 

where T> is as in (JJJ, while 

r =(« p ^{i, 1 5)' (i7) 

and Q = diag(— 1/2, 1/2,0). Although P is a projector, 
(2Q) 2 = diag(l, 1,0), and as a consequence 1Z possesses 
one negative eigenvalue, (1 — y/2)/2, of multiplicity two. 
Any bath for which the Kossakowski coefficient matrix D 
has support only in the negative eigenspace of 1Z would 
generate a negative dt£(0), and therefore entangle the 
initially separated state p(0) = |+)(+| ® |+)(+|- 

A simple explicit example in which this happens is 
given by the following two-parameter matrix T>, with 



/ 1 -ia 
A = C = ia 1 

V 00 












-b 
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where a and b are real constants. |22( Positivity of T>, 
required by the complete positivity of the subsystems 
Markovian dynamics JSJ, is guaranteed by a 2 + b 2 < 1. 
Inside this unit disk, the region for which df £ (0) in (|16|l is 
negative is characterized by the condition a + b > 1. Ac- 
tually, by changing the initial state p(0) and the probe 
vector I ■(/>), one can show that entanglement is created 
in all four disk portions outside the embedded square 
I a ± 6| < 1. Notice that inside this square T> is pos- 
itive, so that there the time evolution of the partially 
transposed density matrix p(t) generated by (|10J) is also 
completely positive: in this case, entanglement can not 
be created for any choice of the initial state p(0) and of 
the vector 

Now that we have shown that a Markovian dynam- 
ics can indeed entangle the two subsystems via a purely 
noisy mechanism, let us discuss in more detail the condi- 
tion for entanglement creation. Although in general the 
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basis vectors |a<), \bi), introduced in are not eigen- 
states of 03, they can always be unitarily rotated to the 
basis |±): 

\ai) = U\+) \a 2 )=U\-), 

\h) = V\+) \b 2 ) = V\-}. (19) 

The unitary transformations U and V induce orthogo- 
nal transformations hi and V, respectively, on the Pauli 
matrices: 



(20) 



With these definitions, for a generic separable initial state 
ijjJJ and arbitrary vector such that £(0) = 0, the con- 
dition dt£(0) < for entanglement formation can be ex- 
pressed as the following expectation value over the prod- 
uct of 6 x 6 matrices: 



w < 



(21) 



where T> is as in fT2")l (with £T< 12 ) set to zero as explained 
before), while the remaining matrices are given by: 



w = 



u 
V 



* = 



V>21 13 







"V>12 13 



(22) 



and the components of the 6- vector w by the Pauli matrix 
elements: 



Wi = (+m\-) 



i = 1,2,3 



(23) 



A more manageable condition for checking entangle- 
ment production can be obtained by noticing that l|21|) 
is quadratic in the components ^12 and "021 of \ip). By 
suitably rearranging the expression in (|21|l , one can then 
show that entanglement is generated if the following in- 
equality, independent from the probe vector lip), holds: 



(u\A\u) (v\C T \v) < \ (u\Ke(B)\v) 



the 3- vectors \u) and \v) are not completely arbitrary: 
they contain the information about the starting factor- 
ized state (JSJ, and their components can be expressed as: 
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Therefore, a given bath will be able to entangle the two 
subsystems evolving with the Markovian dynamics gen- 
erated by © and characterized by the Kossakowski ma- 
trix J7J), if there exists an initial state |ai)(ai| (g> 
or cquivalently orthogonal transformations hi and V, for 
which the inequality (|24|l is satisfied. 

The condition l|24|) can thus be used to check the entan- 
gling power of specific Markovian time evolutions. As an 



example, consider a bath leading to a Kossakowski ma- 
trix (Q) for which A = B = C; this choice corresponds to 
a special case of collective resonance fluorescence 
Provided the hermitian matrix A is not symmetric, one 
can easily prove that there are initial states of the form 
© with \ai) = \bi) that will get entangled by the noisy 
dynamics. Indeed, in this case condition 124(1 reduces to: 



(u\lm(A) \u)\ > 



(26) 



which is clearly satisfied for any \u) outside the null 
eigenspace of Xm(A). When A is real however, (|26|l is 
violated and entanglement is not created, since the par- 
tial transpose state p(t) evolves in time with a completely 
positive dynamics. 



After completion of the manuscript, our attention has 
been drawn to Refs. 0, Hil Eof which have connec- 
tions with the topics discussed in this letter. 
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